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Abstract 

Quantum estimation theory provides optimal observations for var- 
ious estimation problems for unknown parameters in the state of the 
system under investigation. However, the theory has been developed 
under the assumption that every observable is available for experi- 
menters. Here, we generalize the theory to problems in which the 
experimenter can use only locally accessible observables. For such 
problems, we establish a Cramer-Rao type inequality by obtaining an 
explicit form of the Fisher information as a reciprocal lower bound for 
the mean square errors of estimations by locally accessible observables. 
Furthermore, we explore various local quantum estimation problems 
for composite systems, where non-trivial combinatorics is needed for 
obtaining the Fisher information. 



1 Introduction 



In many experimental situations, we are not allowed to have a large num- 
ber of data enough to determine unknown parameters such as coupling con- 
stants of hypothetical interactions. In some cases, the number may be fairly 
small and it is crucial to theoretically explore the best estimator for the pa- 
rameter from the small number of our available data. The problem becomes 
prominent for quantum systems, since optimal estimation must be well recon- 
ciled with inevitable quantum uncertainty arisen from available observables 
and unknown parameters. In such situations, the quantum estimation theory 
can play a significant role; for detailed reviews, we refer to Helstrom pQ and 
Holevoj2]. The theory provides the best observation on the system for the 
estimation with the minimum value of the estimate error. 

Although the ordinary quantum estimation theory is certainly powerful 
for many estimation problems, the theory includes an implicit assumption, 
which is not realistic in some of the practical experiments. The assumption 
is that every observable of the system is available for the observer or the 
experimenter. Contrary to the assumption practically available observables 
are often restricted. For instance, it is a common situation in experiments 
that a particle is contained inside the laboratory at the origin of the time and 
that the experimenter can only use measuring devices inside the laboratory. 
However, according to the time evolution the particle may go out of the lab- 
oratory, so that the ability of estimating the state parameter is restricted to 
measuring devices inside the laboratory for the later time. Another example 
is found in elementary particle physics. It usually happens due to the limit 
of the present technology of measurement that our apparatus can probe only 
low-energy portions of the total Hilbert space with visible signals. Thus the 
observables are certainly restricted. In such situations, the observable pro- 
vided by the ordinary theory for the best estimate may in general not be 
available . Then, the question becomes relevant as to what is the best esti- 
mate among those which are accessible only by use of restricted observables. 
Let us generally call such estimations local quantum estimations. 

In this paper, we elaborate the formulation of quantum estimation theory 
for local quantum estimation problems on an unknown parameter g. For the 
restricted density operators measured by our apparatus, a Fisher information 
is introduced. Then, we prove the quantum Cramer-Rao type inequality 
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for the local quantum estimation for g. The observable is specified that 
attains the equality and yields the best local quantum estimate for g by its 
measurement. It is also pointed out that there exist nontrivial aspects in the 
analysis of the local quantum estimation for the composite system of identical 
subsystems. In that case we have two natural estimations and corresponding 
two Fisher informations for the unknown parameter g. The first alternative 
takes a simple form to apply, but may give a smaller value of the Fisher 
information. The second alternative is able to give a larger value of the 
Fisher information and generates a better estimate for g, but have a pretty 
complicated form to deal with, compared to the first alternative. Especially, 
calculation of the second Fisher information requires solving independently 
evolutions of many descendant operators. 

In section 2, a brief review on the standard quantum estimation theory is 
given. In section 3, we discuss more physically the quantum estimation prob- 
lem, including the biased-estimator case. Several expected advantages of the 
quantum estimation are also reviewed. In section 4, we introduce the notion 
of local quantum estimation problems. In section 5, a quantum Cramer-Rao 
inequality for local quantum estimations is established. In section 6, Fisher 
information is discussed for unnormalized pure states. In section 7, we reveal 
nontrivial aspects of local quantum estimations for composite systems. In 
section 8, two general formulations are proposed for local quantum estima- 
tions for the composite system. In section 9, a formula which is useful for the 
evaluation of one of two sorts of Fisher information for composite systems is 
given. In section 10, we apply our formulation for local quantum estimations 
to a decaying two-level system with a small unknown parameter. In the final 
section, we summarize our results of this paper. 

2 Cramer-Rao bound for quantum estimators 



Let us briefly review quantum estimation theory in this section. For a 
detailed review we refer the reader to Helstrom |T] and Holevoj2]. 

Let S be a closed quantum system described by a Hilbert space H. We 
assume that the Hamiltonian H to t has a constant real parameter g £ G, i.e., 
H to t = Htot(g), where G is the set of possible values of g. The evolution 
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equation for the density operator p tot of S is given by 

ihdtptat = [H t ot, Ptot}- (1) 

Then, the density operator of S at a given time t depends on the time t and 
the parameter g, i.e., 

Ptot = ptot(g,t). (2) 

We shall consider the following quantum estimation problem for the pa- 
rameter g. In order to estimate the parameter g in H tot (g), we assume that 
one measures an observable A at time t, and the output g is taken as the 
estimate for g. Thus, the observable A plays the role of the estimator of this 
statistical estimation problem. 

By the Born statistical formula, the expectation value of the measurement 
output g in the state p ta t(g) = Ptot(g, t) is given by 

E g [A] = Tr[Ap tot (g)}. (3) 

Then, E g [A] is the mean of one's estimate g for the given true value g(E G). 
The variance of the estimate ~g for the true value g(E G) is given by 

V g [A] = (A g Af = Tr{A* Ptot (g)] ~ Tr{Ap tot (g)] 2 , (4) 

where A g A is the uncertainty of observable A in the state Ptot{g)- The mean- 
square error e 2 g [A] of the estimate g for the true value g is defined by 

e 2 g [A]:=E g [(A- g n (5) 

By a simple manipulation, we obtain the relation 

e 2 g [A] = V g [A} + (E g [A]-gf.. (6) 

The estimator A is called unbiased, if the mean estimate is correct for any 
possible values g(E G), i.e., 

E g [A]=g (7) 

for any g(e G). In this section, we shall confine our attention to unbiased 
estimators A. As seen in eqn(jnj), for the unbiased A, the variance V g [A] 
represents the mean-square error of the estimate g for the true value g(& G): 

e 2 9 [A] = V g [A}. (8) 
3 



The lower bound for V g [A] is given by the well-known quantum Cramer- 
Rao inequality as follows. The symmetric logarithmic derivative (SLD) L(g) 
for ptot(g) is defined as a self-adjoint operator satisfying 

dgptot(g) = -j[Ptot{g)L{g) + L(g)p tot (g)]. (9) 

Note here that 

E g [L(g)}=Tr[p tot (g)L(g)} = 0, (10) 

due to the normalization condition of p to t- By the above relation, the SLD 
may not be determined uniquely; however, any two solutions L\(g) and L^g) 
satisfy the the relation |S| 

L\{g)ptot{g) = L 2 (g)p tot (g). (11) 

Thus, the operator L(g)p to t(g) is uniquely determined. The Fisher informa- 
tion J g of the parameter g in ptot(g) is uniquely defined by 

J g = Tr[L 2 (g)p tot (g)]. (12) 

Then, every unbiased estimator A satisfies the quantum Cramer- Rao inequal- 

ityfflU 

4[A] = V g [A] > 1. (13) 
A simple proof is given in the appendix. 
3 More Physical Review of Quantum Estimation 



In this section, we shall discuss real experimental procedures for esti- 
mating the unknown parameter g in the framework of quantum estimation 
theory given above. In real experiments, rigorous unbiased estimators that 
satisfy eqn(J7j) globally in the parameter space of g are usually not available 
for experimenters from technical reasons. Instead, only biased estimators 
are available. Even in such real situations, as we shall show in the following, 
the quantum estimation theory described in the previous section plays an 
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active role. We also discuss some advantages of quantum estimation theory 
to provide more efficient methods in several estimation problems in physics. 

We shall first consider a quantum estimation for a parameter g in H to t(g) 
performed by a measurement of a general observable A for a single sample 
at time t. Here A is not assumed unbiased. Assume that we get an outcome 
a in the measurement and we make the estimate g for g as a function of the 
outcome a. Usually, this function g = 0(a) is obtained by the following way. 
The relation 

E g [A\ = a (14) 

between the true value g and the mean output a from the measurement can 
be often solved theoretically as a function 

9 = 0(a)- (15) 

Usually, sensible experiments are designed to possess a suitable domain C?(C 
G), which includes interesting values of g, and to allow the relation 

2 = 0(a), (16) 

applying the above function to the measurement output a, gives a good 
estimate g for the given true value g(G Q) from the output of the single 
A measurement. Note that, due to the quantum nature, the observable A 
generally possesses nonzero values of the variance V^L4]. The variance has 
a close relationship to the problem; to what extent the estimate g can be 
trusted. For example, if V g [A] is extremely small, then, even in the single 
measurement, the observed value a must be almost equal to the expectation 
value -5^1/4] for the correct value of g. Consequently, the estimate g has 
to almost coincide with the correct g value. From the viewpoint of the 
dimensional analysis, it is rather straightforward to introduce an expected 
error 5<7 ff [A] of the estimate ~g~g[A] for the true value of g as follows. 

Now, we shall give a justification of the above estimate and error eval- 
uation from the viewpoint of quantum estimation theory. When V g [A] is 
small enough, a domain Q in the space G of possible values can be chosen so 
narrow that E g [A] is linearly expanded around a physically interesting value 

g (egy. 

E g [A] = E 9o [A] + d g E g [A]\ g=go (g -g ) + ((g - g f) . (18) 
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Here it is quite useful to remind that a single measurement of the observable 
A simultaneously implies a single measurement of an observable f(A), where 
f(x) is an arbitrary real function of x. If an output a of the observable A is 
obtained, it is interpreted that an output /(a) is observed for the observable 
f(A) in the same measurement. In what follows, in order to make a useful 
choice of /, let us impose on the function / the locally unbiased condition: 

E g [f(A)]=g + 0((g-g ) 2 ). (19) 

Then it is noticed that eqn (|19|) is satisfied for g 6 Q by a linearized function 
f(x) such that 

f(x) = e +q - Mm (20) 

A ] d g E g [f{A)\\ g=g J 90 d g E g [J{A)\\ g=g ; 1 > 

By linearity of /, for the mean output a = E g [A] we have 

f{a) = f(E g [A]) = E g [f(A)), (21) 
so that from eqn ()15|) we have 

( J ) (a) = f(a) + 0((g-g ) 2 ). (22) 

Thus, g = 0(a) is now reproduced by substituting the output a into the 
function f(x). By use of eqntJHJ and eqn()19p. the mean-square error of the 
estimate g as the output /(a) of the f(A) measurement is evaluated as 

t 2 9 [f{A)} = V g [f{A)]+0({g-g y). (23) 

From eqn(j20J), the variance of f(A) is evaluated as 

V g [f(A)} = p V ? [A \ A „ 2 = {5g g [A]f (l + 0(g- 9o )) , (24) 



thus, it is verified due to eqn (j23j) and eqn (j2Tj) that 

e 2 9 [f{A)] = {5g g [A])\l + 0{g-g )). (25) 

Hence the validity of the error evaluation by (og^A]) 2 definition has been 
shown. 
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It is a significant result from quantum estimation theory that even for the 
biased observable A, the quantum Cramer-Rao inequality can be proven: 

Here the Fisher information is defined by eqn()12|) and the equality can be 
achieved by taking A oc L(g) for each value of g. This can be shown by 
adopting not A but the local unbiased operator f(A), and returning to the 
general argument for the unbiased case in section 2. It is also possible to 
prove by using the biased observable A straightforwardly. The proof can be 
seen in the appendix. 

Next we shall discuss physically relevant cases, the N-samples systems, 
by naturally extending the single-sample argument. Let us take a composite 
system which consist of N identical S subsystems. Assume here that the 
density operator is independent and identically distributed; 

p£ t ) (9,t) = p(g,tf N . (27) 

Now the average estimators defined by 

1 N 

A(N) '■= ^E 1 ®'--® 1 (28) 

are available. According to the quantum law of large numbers the mea- 
surement data for A^ are going to be normally distributed with the average 

1 /2 

-E^ [A] and the standard deviation ( Va j^ ) when the number N becomes 
large. Since the expected error is solely a pull back of the quantum deviation 
of the observable, we can trust the estimate, 

g = g±6g B [AW], (29) 

in the 1-cr precision for the large number cases. 

Now let us discuss the Fisher information for the N-samples cases. The 
SLD for the composite system ZA^ are defined for general density operators 
p\ot by 

d g pS = l[pSLW+LWp£ t \ (30) 

For independent and identically distributed (i.i.d.) density operators, it is 
easily derived that the composite SLD L^ N ' is given by 

N 

LW = ^1®.®1®I®1®.--®1, (31) 
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where 1 is the identity operation and L is the SLD for the subsystem defined 
by eqn®. This result yields the following simple relation for the Fisher 
information of the composite system = Tr[(L^ 7V ^ ) ) 2 p^]: 

jW = NJ V). (32) 

By virtue of the Cramer-Rao inequality, it is easily noticed that the op- 
timized estimation for g in a single measurement of the composite system 
is achieved by adopting the average estimator iP^ = j^L^ N \ The expected 
error is given by 

Sg 9 [L {N) ] = ~r±==- (33) 

This coincides with the usual error of the estimation for g based upon N 
independent data of measurements of L^ 1 ' for the N subsystems. However, 
stress that there is no need to measure N times the estimator iS 1 ' for each 
subsystem S to achieve the estimate. In the quantum estimation, just one 
measurement of the single observable L^ N ' yields the best estimate. Other 
relative- difference components like L ® 1 • - • Cg) 1 — 1 ® L •-•(£) 1 remain un- 
measured. This saves effectively the number of processes in the estimation 
and exposes an advantage of the quantum estimation. 

When the entanglement between the subsystems is available, it is possible 
[2] that the large N behavior of 5g g [L^] can be improved beyond the 1/y/N 
factor as 

6g g [LW] cx 1 (34) 
This reveals another advantage of the quantum estimation. 

4 Quantum Estimation by Local Observables 



We shall now consider the following constraints on the quantum estima- 
tion problem discussed above. In the above general formulation, we have 
assumed that every observable A of the system S is available for our mea- 
surement to fix the g value. However, in practice the available observables 
are restricted. For instance, it is a common situation in experiments that a 
particle described as the system S is contained inside the laboratory at the 
origin of the time and that we can only use measuring devices inside our 
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laboratory. However, according to the time evolution the particle may go 
out of the experimental apparatus or our laboratory, so that for the general 
t, our ability of estimating the parameter g is restricted by the measuring 
devices inside the laboratory. 

Let M. be a subspace of 7i. The projection of H onto M. is denoted by 
P. In this paper, we consider the following two constraints. 

(i) The initial state is supposed to be supported by P, i.e., 

p tot {t = 0) = Pp tot {t = 0)P. (35) 

(ii) The available observables O for our measurements are restricted to 
those of the form 

O = PXP + y(I - P), (36) 

where X is an arbitrary observable on 7i and y is an arbitrary real number. 

Let {\a)} be the orthonormal basis of M. and {|«)} be the orthonormal 
basis of 7i extending {|a)}, i.e, {\a)} C {|ct)}. Then, condition Eq. (|33j) is 
equivalent to the relation 

Ptot(t = 0) = £ \a)(a\ Ptot (t = 0)\a')(a'\. (37) 

a,a' 

Thus, the density operator initially has only matrix elements inside Ai, and 
according to the time evolution, the density operator p tot may have matrix 
elements outside of M.. 

In the case of estimating the parameter g by observing a particle S initially 
localized in a box using the measuring devices effective only inside the box, 
the subspace M. corresponds to the space of wave functions localized in the 
box. In this case, the assumption that the particle is initially localized inside 
the box is represented by condition (i). Since we assume that we know that 
the particle inside the box at the origin of the time, by measuring later, for 
instance, the weight of the box , we can measure I — P via a negative result 
of the measurement. Since the measuring devices are only effective inside the 
box, the measuring interaction couples only with the observable of the form 
PXP so that it is natural to assume that they can measure only observables 
of the form PXP. Therefore, the set of available observables are considered 
to be restricted to those given by condition (ii). 

Initially the density operator p tot have only matrix elements in 
However, in the course of the time evolution, p tot can have matrix elements 
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outside of {|a)}. For any nonnegative time t > 0, we define the accessible 
density operator p\\ (t) for the subspace M. by 

p\\(t)=Pp tot (t)P- (38) 
Obviously, p\\ has the matrix representation 

p\\(t) = [(a\p tot (t)\b)] ■ (39) 

Then, by the corresponding properties of ptot(t), the operator p\\(t) is positive 
and satisfies 

< Trp|| (t) < 1. (40) 

In what follows, we shall consider the time domain of t from t = to the 
time just before t = t* such that Trpy(t = t*) = 0, where we allows the case 
= oo, so that we have 

0<Tr P ||(t)<l, (41) 

for t G [0,Q. 

From condition (ii), the available estimators A on Ai are naturally re- 
stricted and satisfy the relation 

A = A|| + o_l(1- P) = A\ (42) 

where 

A\\ = PA\\P. (43) 

Using the definitions in eqns (|38|h (j4^|) and it can be shown that the 
expectation value of the available estimator A is given by 

(A) = Tr( Ptot (t)A) = Tr( P ||(i)A||) + a ± (l - Ttp„(*)). (44) 

In order to define rigorously the notion of the "local" estimators A corre- 
sponding to the available estimators A in the restricted situation, let us ex- 
tend the space M. to a one-dimension-larger Hilbert space M. by adding to the 
basis of M. a normal vector \B) orthogonal to every |a),i.e., A4 := A4®C\B) . 
Then, the local estimators A acting on A4, which corresponds to the available 
estimator in eqn()42)l. are defined by 

A = A\\ +a ± \B)(B\. (45) 
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In particular, note that 



1-P=\B)(B\. (46) 

Since the state \B) represents the inaccessible states by our local observation 
as seen in eqn (|4fij) . we call \B) the blank state. 

Further let us introduce the local density operator p acting on M. and 
corresponding to p\\{= PptotP) by 

p = p\\ + (l-Tr P]l )\B)(B\. (47) 

It is easily seen that p is positive and of unit trace. By a simple manipu- 
lation, the expectation value of the available estimator A in eqn (|44|) can be 
reexpressed by use of the local estimator A and the local density operator p 

as 

(A) = Tr(pi). (48) 
5 Cramer-Rao Bound for Local Quantum Estimators 



In what follows, we shall consider the quantum Cramer- Rao inequality for 
the quantum estimation problem for the coupling constant g in the Hamil- 
tonian H tat (g) by using only local measuring devices. By the time evolution, 
the local density operator p = p(t, g) introduced in the previous section de- 
pends on the time t and the parameter g. Now we assume that one measures 
a local estimator A at time t, and the output a determines the estimate g~a[A] 
for g via the relation: 

E MA] [A}=Tr[Ap(t,g- a [A])}=d. 

Stress that, in this estimating process of g, we are allowed to use only local 
estimators A instead of arbitrary observables in the theory. 

The variance of the local observable A for the correct g value is certainly 
given by 

V g [A] = Tr[A 2 p(t,g)] - (Tr[Ap(t,g)}) 2 . (49) 
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Then it is required at a given time t to find the minimum value of the expected 
error defined by 



Sg a [A] 



1 J (50) 



\ (d g E g [A])* 



It is shown that this problem is resolved by use of a solution of the problem 
on the estimate for the parameter g by arbitrary observables O on Ai as 
follows. 

We define the local SLD L(g) on M. for an arbitrary local density operator 
p(g) in eqn p7|) self-adjoint operator satisfying 

d g p(g) = ~(L(g)p + pL(g)), (51) 

L(g) ] = L(g). (52) 
Since Trp(g) = 1 for any g G G, we have 

Tr[L(g)p(g)} = 0. (53) 

It is easy to construct a solution of eqn (Jf)T|) by introducing a SLD operator 
L(g) on M. for the accessible density operator p\\. The SLD L(g) on M. is 
defined by 

dgP\\ = ^{L{g)p\\+P\\L{g)), (54) 

L(gy = L{g), (55) 
PL(g)P = L(g). (56) 

Due to the fact that Pp\\P = p\\, we can find, at least, a solution of eqn (|54j) 
for the SLD with PL(g)P = L(g). Once the SLD L(g) is given, then it is 
proven by a simple algebra that the operator defined by 

L(g) = L(g) + <9,ln[l - IVp||G/)]|fl)<fl| (57) 

satisfies eqns (Jo"Tj) and ()52p. thus it is a SLD on M. for p(g). 
Here it is carefully noted that we may have 

Tr[L(^)p||(^)]^0, (58) 

since the trace of p» (g) is not necessarily normalized. 
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The operator L(g) is determined uniquely up to the support of p(g); 
any two solutions Li(g),L 2 (g) satisfy L 1 (g)p(g) = L 2 (g)p(g). The Fisher 
information J g of the parameter g in p(g) is uniquely defined by 

[Tr[L(o)pi|(o)ll 2 

J g = ^[L{gfp{g)\ = Tr[L(g) 2 P]l (g)} + L , l? 1 )' , (59) 



TrhG?)] 

where we have used eqnslHTJ) and (fHTjl . Then, for the arbitrary observables 
O on Ai we have the quantum Cramer-Rao inequality 

where J g is given by eqn ()59|) . In order to apply the result to our local es- 
timator problem, it is crucial to notice that the SLD in eqn ()57|) takes the 
precise form of the local estimator on Ai in eqn (J45|) . Therefore, the equal- 
ity can be attained by a local estimator. This indicates that the following 
quantum Cramer-Rao inequality for arbitrary local estimators A on Ai re- 
ally holds for the local density operators p corresponding to the accessible 
density operators p\\(g)(= Pp to t{g)P)- 

where the Fisher information J g is given by eqn (JoU|) . For a given g G G, the 
equality is attained by a local estimator A Q (g) such that 

A M oc Ut) = m - T ^mm. [B){ B\. (62) 

Note that the local estimators which give the minimum expected error such 
as A Q are unique only up to a factor and an additive term proportional to 
the identity operator. For instance, an estimator such that 

A M oc L (9) + ™, (63) 

which has no matrix element for the blank state, also attains the equality. 
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6 The Fisher Information for Unnormalized Pure States 



In physics, it often happens that the measurement device is able to probe 
only a small part of the physical states of the total system. Even in such sit- 
uations, non-unitary formulations are sometimes available. The state vectors 
\^{t)) are governed by equations of motion with non-Hermitian Hamiltoni- 
ans and evolve deterministically in the subspace A4, which is accessible by 
the experimental devices. Such examples are found in the various fields of 
physics, including the scattering problems with weak absorption of quanta in 
the nuclear physics and the quantum optics, the flavor-oscillation studies in 
the elementary particle physics and so on. The information about the cou- 
pling constant g in the equations of motion is imprinted on the state vectors 
\ty(t,g)) during the time evolution. 

Let us evaluate the Fisher information for the pure state \^>(t,g)). The 
accessible density operator for the pure state reads 

p {l (t,g) = \#(t,g))(#(t,g)\, (64) 

where Tr[p||(0, g)] = 1 and at an advanced time £(> 0) the following relation 
holds: 

0< < 1. (65) 

We define the SLD operator L on A4, in the same way discussed in the previ- 
ous section, for the accessible density operator py. Note that the operator L 
is not uniquely determined due to the purity of py , however, the ambiguity is 
not relevant at all for the Fisher information, as commented in the previous 
section. It is shown that we have a SLD, 

<? Tr[0 gP ||] 

Ap\\~ mvrf i \2 P\h ( 66 ) 



Tr[p,|] ^" (Tr[p||]); 

as a simple representative and the Fisher information itself is uniquely eval- 
uated by 

j = 4 ( o w *> - MM?) + 4 IJfe(*IWI 2 (67) 

where \d g ^) := d g \^f(t, g)). This result is an extension of that in the reference 
[3], where the normalized pure state theory is analyzed. The relation enables 
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us to evaluate easily the Fisher information for many unnormalized pure 
state theories. 

In eqn (167)1 . one may worry about the apparent divergence of the third 
term at (^l^) = 1, because the state evolves initially from the normalized 
state. However, for ordinary physical systems, the early behavior of the norm 
(^l^) is given by 

(*(t,g)\y(t,g))~l- a (g)t 2 , (68) 

where a is a positive function of g. Thus the third term is evaluated in the 
early era as 

1 -<*|*> a(g) *■ [bJ) 

Hence, the limit t — ► +0 of eqn ()67j) exists without any problems. 

7 Problems of the Composite System 



In the local estimation problem, some nontrivial aspects appear in the 
composite system analysis. Suppose a system S. Let us assume our mea- 
suring device for S is able to access only a subspace M. of the Hilbert space 
of S. Later let P denote the projection operator onto A4, &ndD M denote 
the dimension of the subspace M.. The accessible density operator on M. 
is denoted by p\\. The operator p\\(t, g) evolves in the subspace M. and be- 
comes dependent on the coupling constant g in the equation of motion. Let 
us consider a composite system S® N composed of N identical S subsystems. 
For instance, suppose that an independent and identically distributed (i.i.d.) 
initial condition is set for the total density operator pt o t{0) of the composite 
system. Also assume that the unitary evolution of the total system is fac- 
torized, i.e., U^ N \t) = U(t)® N . Even in such a simple situation, it can be 
pointed out that we have, at least, two natural alternatives for the estimation 
of g as follows. 

The first alternative is rather simple. In the procedure, one firstly cal- 
culates the accessible density operator /?f| = P® N p tot P® N for S® N which is 
reduced to a direct product defined by 

pW(t,g):=pf N (t,g). (70) 
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A local density operator for the accessible density operator p^ can be de- 
fined straightforwardly by 

p[ N) :=p\P + (l-Tr P \P)\B)(Bl (71) 

where \B) is the blank state. Let denote the Fisher information based 
upon the first local density operator p^\ 

The estimation problem in the composite systems is nontrivial because 
we may have a construction of another local density operator for S® N . We 
are able to define at first the local density operator for each subsystem S. 
For the i-th subsystem Si, the local density operator pi corresponding to p\\i 
is written as 

Pi = PIN + (1 - Trp |K )|^>{^|, (72) 

where \Bi) is the blank vector for the i-th subsystem S«. Then we can define 
naturally the second local density operator p^ for the composite system 
S 07V by a direct product as follows: 

pf } :=rW (73) 
i=i 

Let denote the Fisher information based upon By construction, 

the Fisher information j( N ) for the i.i.d. density operator is calculated as 

j(A0 =A rj(i). (74) 

As seen above, there exist two independent Fisher informations for the 
composite system. Then, it is an important question; which alternative of 
the formulation gives us a more precise estimate for g, that is, which Fisher 
information, j( N ) or j( N \ is larger than another. The problem should be 
addressed for the general initial conditions for the density operators, beyond 
the above i.i.d. situations. Note first that the operators p[ N ^ act on the 
Hilbert space M m © C\B) and the dimension of M® N © C\B) is given 
by D 1 = (Dm) N + 1- On the other hand, the operators p^ act on the 
Hilbert space Ai® N and the dimension of Ai® N is given by D2 = (Dm + 1)^. 
Since D 2 > D\ always holds, it is naively expected that the second Fisher 
information j( N ) is not less than the first Fisher information j( N \ This 
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guess can be proved affirmatively by use of the monotonicity argument for 
the Fisher information as will be mentioned later. 

The above argument has been limited to the i.i.d. cases. In order to an- 
alyze the composite-system estimation generally, we must extend the above 
two formulations. Especially, nontrivial analyses are required to define J^ N \ 
These are formulated in section 8. 

In section 9, it is also pointed out that evaluation of the larger Fisher 
information J^- 1 requires solving time evolutions of various density operators 
corresponding to different initial conditions. Such a feature does not appear 
in the evaluation of both the standard Fisher information in the usual cases 
and the smaller Fisher information the local estimation. 

8 General Formulation for the Composite System 



The available estimators for the composite system S® now reads 

= Yl u kl ... kN A kl ® • • • ® A kN , (75) 

which is just a natural extension of eqn([42|) . Here A k denote the available 
estimators for the subsystem S, which take the form in eqn (}4"2|) and ui kl ... kN 
are real coefficients. The corresponding extension of eqn (J45|) is also possible. 
The local estimator A^ N ' corresponding to the available estimator A^ N > is 
defined by 

iW= u) kl ... kN A kl ®---®A kN , (76) 

where A k are the local estimators corresponding to A in eqn (|45l) . 

In order to define the two Fisher informations and beyond 
the i.i.d. condition, let us consider the most general local density opera- 
tor p£?(0) = P m p { t o t \0)P m as the initial total density operator. In the 
unitary time evolution of the total system, 

p£ t \t, g ) = uW(t, g ) P £ t \0)uW\t,g), (77) 

the total density operator becomes to have matrix elements between the 
inaccessible states. 
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Even for the general initial conditions, the definition of the first Fisher 
information jW is essentially unchanged. Let us introduce the accessible 
operators p\\ by reducing the total density operator p£2 as 



S.N) _ p ®N AN) p <g)N 
'\\ — r Ptot r 



(78) 



For the accessible density operator pi , a SLD operator L,W is defined by 



(79) 



( L (W))t = l (jv) j 



(80) 
(81) 



According to eqn(|59|l. the Fisher information 

j(N) ig 

defined straightforwardly 

as follows. 



Tr 



(N) 



(Tr 


L (N) 


^11 


r 


1 - 


-Tr 


P \\ 





(82) 



Next, in order to define the second Fisher information j( N \ what we want 
is a proper definition of a local density operator p( N ' acting on Ai® N such 
that the total density operator p^ is reduced into p^ N \ Here, it is quite 
natural to impose that expectation values of all the available observables 
AW for p\ot are equivalent to those of the corresponding local observables 
AW for P W : 

Tr[AWp tot ] = Tr[AWp^}. (83) 

By some manipulations it is soon noticed that the above constraint is really 
satisfied by defining the local density operator 

pW as 

follows. Let index a,- 

for j = 1, . . . , N below take index aj for states in Ai or the index B for the 
blank state. Then the matrix elements of pW on AA® N are given by 



(0L1OC2 • • • ttjvlP' 
JV 



IV 



x e n [^B^ fe i + (i 

Xl—X}f k=l 

xTr[(P 1>Xl ® P 2 ® ■ • • ® 



- <J« fc fl)<J<r fc o] II [<^,lA V l + 
fc'=l 

-P/V,xjv)Ptot]; 



(1 - 



54) 



18 



where Tr stands for the trace operation on 7i® n . For m = 1,...,N, the 
subscript x m takes or 1 and the operator P m ,x m is defined by 

P m ,x m = \a' m )(a m \, if x m = (85) 
P m>Xm = l-P, ifx m = l. (86) 

By construction the Hermicity of the operator pW is trivial. Further, taking 
A W = 19N 

in eqn (|83|) yields the normalization condition: 

Tr[pW] = 1. (87) 

The positivity of pW is also proven as follows. Suppose an arbitrary vector 
on M m : 

1^) = E C ai ... aN \ai - ■ -a N ) 

= E \*Hvi h] ), (88) 
(u— »*) 

where 

l*[0]) = E C'ai-ajv K •■ -ajv), (89) 

l*[l]>= E CBaa-awlSaa-'-ajv) (90) 

and so on. Then, using the definition of p^ in eqn (J84)) . the expectation 
values of p( N > for the arbitrary state vectors |\&) are evaluated as follows. 

= E Tr[iV, fe]P W] 

= E Tr[P [n .., fc]Pio4 ], (91) 

where P^...^] = \^[i v -.i k ])(^[i v -.i k }\ and P^...^] are defined by replacing \ B)(B\ , s 

in the operator P[i v ..i k ] to 1 — P. Noting that the operators P^...^] can be 

expressed as P[i r .. ifc ] = E/3 l$/J,[»i-i fc ])($,9,[ii-t k ]l b Y use of vectors 1$^...^]) 
in the total Hilbert space, it is proven that 

<*|P (A °I*>= E (^ 1 [i 1 -^]bS ) |^ 1 [i 1 -y>>0. (92) 

P,h,—ik 
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Taking account of the normalization condition in eqn (J87|) . this implies the 
positivity of the operator p( N \ 

Since Tr[p^] = 1, we can define in the usual way a SLD operator C for 
the local density operator p( N >: 



9gP 



(N) 



Then, the Fisher information is defined by 

J(N) = Tr UN) C 2 



(93) 



(94) 



Now let us comment on the inequality > j ( N > , using the monotonicity 
of the Fisher information. The point is that there exists a mapping R of 
the density operators defined in A4® N onto the density operators defined in 
Ai® N © C\B). Let us denote P\\ a projection operator onto the subspace of 
vectors in J\/[® N that do not include the blank states at all. Denote P± a 
projection operator onto the subspace of vectors that include more than one 
sub-blanck vectors \Bi). It should be noted that 



P\\P {N) P\\ = P\\ 
Let us define the mapping R as follows. 



(AT) 



(95) 



R{p 



P\\P {N) P\\ +Tr \P±p {N) ] \B)(B 



P\\ 



(N) 



+ Tr 



P ±P W] \B)(B\. 



By definition, it is clear that the mapping is linear and of unit trace: 



Tr 



R[p 



(N)] 



Tr[P|| P W]+Tr[P ± p 



1. 



(96) 
(97) 



It is also easily seen that this mapping is completely positive, since so are 
p (N) ^ P IIP ( N )P {1 and P W Tr [P ± pW . 

Using the relation Tr[p^] = 1, we obtain 



R[pW] 



Tr\pW])\B)(B\. 



(98) 



Then, the first Fisher information j( N ^ is given by Tr[^^)(L) 2 ], where L is 
the SLD operator corresponding to . According to the monotonicity the- 
orem for the Fisher information jB], it must be satisfied under the projective 
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mapping R that j( N ) > j^ N \ This result does not depend on whether the 
total density operators p^ are factorized or entangled. 

It is worth noting that the information j( N ) possesses a decomposition 
representation. Let us consider an arbitrary subsequence (ii,i 2 , • • • , in) of the 
sequence (1, 2, 3, • • • , N). Define that p^, ...g is a (dimM) N ~ n x(dimM) N ~ n 
matrix which is composed of components of p^ with a im = a' im = B for 
m — 1, • • • , n. The followings are examples. 

(a 2 a 3 ■ ■■a N \p [1] \a! 2 a'z---a' N ) := (Ba 2 a 3 ■ ■ ■ a N \p {N) \Ba' 2 a' 3 ■ ■ -a' N ), (99) 

(a 2 a 4 ■ ■ ■ ajv|p[i,3] |«2 a 4 ' ' ' a jv) := (Ba 2 Ba 4 ■ ■ ■ aN\p^\Ba 2 Ba' 4 ■ ■ ■ a' N ), 

(100) 

(ai|p [2 ,3,.,iv]|a / 1 ) := ( ai BB ■ ■ ■ B\p^\a[BB ■ ■ ■ B) . (101) 

Note that the empty subsequence corresponds to the accessible density 
operator: 

(aia 2 ■ ■ ■ ajv|p[0]|a'iO2 • • • a' N ) = {a\a 2 ■ ■ ■ o,n\P\^ \a[a 2 ■ ■ -a' N ). (102) 

By definitions the local density operators p^,^,-,^] are non-negative, i.e., 
P[i u i 2 ,-,i n ] > 0. For each P[i u i 2 ,-,i n ], we can introduce a partial SLD operator 

£[ii,i2,-,i n ] aS 

dgP[ii,i 2 ,— ,in] = 2 (^[h^,— , in] P[h,i2,—, in] + P[h,i2, — ,in]£ , [h,i2, —,»«])■ (103) 

Then it is possible to rewrite the second information as 

J"" = £ •©,*,. (1«4) 

(il,l2, •■■,ln) 

^[il,i,-,»n] = ^ [^[il.*2.-.»n]^fil,i2,-,»n]. ' ( 105 ) 

Here S(u,i2,-,i n ) means the sum over all the subsequences i 2 , • • ■ , i n ) °f 
(1, 2, 3, • • • , JV), including the empty subsequence 0. The decomposition rep- 
resentation makes the evaluation of easier in many practical applications 
by using a useful formula for the operators P[i u i 2 ,-,i n ] in the next section. 

9 Evaluation of the Local Density Operator 
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The accessible density operators pj| ' can be followed by our apparatus, 
since the operators p^ are completely local by definition. Meanwhile, the 
local density operator p( N ^ has been so far defined based upon the total 
density operator p^P in the previous section. We must say that the definition 
is too formal from the the practical viewpoint, because we seldom know global 
information about the total density p^P due to the limitation of our ability 
to measure the system. For the realistic evaluation of j( N \ it is convenient 
to write down p^ explicitly in terms of locally accessible quantities just 
as the operator p\^ ■ Such a reformulation can be realized for the cases 
with factorized evolutions, i.e., U^ N \t) = [U(t)] (g,N as follows. It should be 
emphasized that we do not need to assume the i.i.d. condition for the initial 
density operator. 

Suppose that a composite system S® N of N identical S subsystems is 
governed by a unitary evolution and that the evolution is factorized for each 
subsystem S: 

p£V) = MtT N p£ t \o)[u(t)r N , (106) 

where U(t) is the unitary time evolution operator for S and pt o t(0) is ar- 
bitrary initial density operator, which may have entanglement between the 
subsystems. 

Let O m := {e a \el = e a , Pe a P = e a } denote the complete basis of the 
available observables acting on M for each subsystem S. Even in our local 
experiments, we are able to define and measure the projective evolutions for 
the available observables e a . The projective evolutions are given by stochastic 
mappings T(g,t)[e a ] which are defined by 

r(g,t)[e a ] :=PU(t)e a U(tfP. (107) 

In various physical systems, the dynamics is first given by not T(g,t) but a 
Lindblad differential form given by 

dtP\\ = T g [ Pll ] (108) 

for arbitrary density operators p\\ on M.. Here T g is a time-independent 
Lindblad super-operator. Then the super-operator T g is related formally to 
the stochastic mapping T(g,t) via 

T g = d t T(g,t = 0). (109) 
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By integrating formally as T(g, t) = e iTfl , we can recover the stochastic super- 
operators T(g,t). 

Stress that the operators T(g,t)[e a ] are completely local quantities we 
can observe. Moreover, the projective evolutions for the composite available 
observables are also completely local quantities, which are written as 

T(g, t)® k {e ai ® • • • ® e ak ] = P® k U{t) m {e ai ® ■ • • ® e ak ){U(t)^P m . (110) 

Our aim in this section is to express the operators P[i lt i 2 , ■■■,{„] of the local den- 
sity operator p( N > only in terms of the accessible operators like T(g, t)® k [e ai ® 
•••®e J. 

Since the initial density operator satisfies p^t (0) = P^p[^(0)P^ N \ the 
operator pioti®) can be expanded using the basis {e a }: 

Ptot(0) = EE' , 'E C «»-«/«i ® e a 2 ® ■■•8>e 0jv , (111) 

ai a2 ajv 

where the real coefficients C aia2 ... ajv is uniquely determined by p£t (0). Af- 
ter rather straightforward calculations, we argue that the following relations 
really hold: 

n 

m=0 (il,-iJm)C(ii,-,i„) 

Xr &(il,i2,-,i»)-(il,-,im) K^X]^'*) [^(il.-JmjbS^O)]] • 

(112) 

Here J2(j 1 ,---,j m )c(i 1 ,---,i n ) means the sum over all the subsequences (ji, • • • ,j m ) 
of the sequence (zi, • • • ,i n ). Tr^ 1) ...j m ^ is a trace operation in terms of the 
(ji; ' ' ' ,3m) degrees of freedom. Tr^ u i 2t ... tin j^(j l! ... t j rn ) means a trace operation 
in terms of the complementary subsequence to the subsequence (ji, ■ • • ,j m ) 
of (ii,---, i n ). When m — 0, Tr(j lr . t j m ) is reduced into the identical opera- 
tion. The operation r^^""] is the time-evolution operator for all degrees of 
freedom removing the (ji, ■ ■ ■ ,j m ) part. Here it is better to note again that 
even though the formula include subtractions, all the operators P[i lt i 2 ,—,i n ] are 
non-negati ve gfr ,^,...,^] > 0. 

In eqn (|112|) . note that the operators P[n,i 2 ,-,jn](0 cannot be evaluated 
only from the knowledge about the accessible density operator p\^\t). It is 
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required to solve independently time evolutions of many descendant opera- 
tors, 



(N-m), 



(g,t) Tr Ul ,., jm) [pS(0)) 



\jl,-,3m] 1 

Here let us just draw the outline of the proof using the case. Substi- 
tuting eqns(|nnj) and (jHU), the following manipulation is possible: 



(a 2 a 3 ■ ■ ■ a N \p {1] \a' 2 a' 3 ■ ■ ■ a' N ) 
(Ba 2 a 3 ■ ■ ■ a N \p( N ^\Ba' 2 a' 3 ■ ■ ■ a' N ) 



= (a 2 a 3 ■ ■ ■ [(1 - P) ® P^ N ^pS(l - P) ® P®^ 

Moreover we can rewrite the equation as follows: 



\a' 2 a' 3 



■a' N ). 



P[i] = Tr i 



P^ N -Vp£ t \t)l®P®( N -V}-Tr 1 \r® N (g,t)[pS(0)} 



It is noticed that the first term of the r.h.s. is reduced using the expansion 
of 

PtoiP(O) as follows. 

= E C aia2 ... aN Tr 1 [U(t)e ai U\t)](PU(t)e a2 U\t)P) ® • • • ® (PU(t)e aN U\t)P) 



^ ] Caia2---ajyTri [e ai ]e a2 



Consequently we arrive at the relation for pm in eqn ()112j) : 



(113) 



P[i]=r^ iV - 1 )(^t)[Tr 1 [p 



(TV) 
tot 



Tri 



(114) 



The proofs for the other components in eqn ()112|) can be achieved in the 
similar ways. 

The relation in eqn ()112|) makes the evaluation of possible, only based 
upon our local knowledge. 



10 A Decaying Two-Level Model with a Small Unknown Param- 
eter 
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In order to demonstrate our formulation explicitly, let us consider a sys- 
tem including a small unknown parameter g. In many physical systems, 
the estimation of such a small parameter often provides significant physical 
information. For example, tiny coupling constants in the elementary par- 
ticle interactions produce only quite rare processes, however, the analyses 
give a lot of important constraints of high energy features beyond the to- 
day's accelerator technology. For simplicity suppose a decaying two-level 
system including the small coupling g. The model has been frequently used 
in physics, for instance, to analyze the flavor-oscillating phenomena in the 
Kq — Kq system j.7j. The Hamiltonian of the example is given as 



H = -ih 



r+ o 
o r_ 



+ gh 



1 

1 



(115) 



where r± > 0, T + ^ T_ and \g\ <C T±, \T + — In the two-level subspace, 
time evolution of the density matrix p» is governed by the following equation 
of motion: 

ihdtp\\ = H P\\ ~ P\\ R] ■ ( 116 ) 

Define the states |±) as 



<7 2 |±) = ±|±). 



117) 



Here <j z is the z-component of the Pauli matrix. Let us estimate the time 
evolution of |±) in the first order of g. It is solved as 



\±(g,t)) = e- r ± t \±)+ i gd(t)\ T ) + 0(g 2 ), 
where the functions d(t) is given as 



d(t) 



e -r+* _ e -v-t 



r+-r_ 

Then the information J± for the state | ± (t)) is evaluated as 

J±(t)=4d(t) 2 + 0(g). 
The Cramer-Rao bound is always achieved by an observable 

A{U) =a y + 0(g), 



118) 



119) 



(120) 



121) 
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where a y is the y-component of the Pauli matrix. 

In this simple model we are able to optimize the measurement time t. 
The Fisher information takes its maximum value in the lowest order 




Jma X = ^ jhj - Ur 1 \ +0(9), (122) 



at 

lnr + — lnT_ ^. . . 
t*= r + _ r + 0(<?)- (123) 

In eqn f)12H|) . when T_ is much smaller than r + : 

r_<r+, (124) 

i r 

£* ~ — p-j- In p^- becomes larger logarithmically. In fact, the time t* can be 
late until the first order estimation breaks down, the time t a ~ — ^-ln^. 
Interestingly, at t = £*, the survival probability for the state | + (g, t)) in the 
two-level subspace is estimated as 

+ (»,*.)> ~ (j^) «1- (125) 

Against a naive expectation, this indicates that the best quantum estimation 
can be achieved at the time after the state has almost escaped from the two- 
level subspace. 

Analysis of a composite system of the two subsystems may be also in- 
structive. Let us first take the initial state as an i.i.d. state, 

l*(0)> = ! + +)• (126) 

For the density matrix p| (t) = \ + +(t))(+ + the fisher information 
defined by 

,*> := ^(LMfi 4- ^^'J (127) 

1 — lrpy ; 

is evaluated as 

j( 2 )(t) = 8d(t) 2 e~ 2r+t + 0(g). (128) 
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Because j^ 1 ' = 4c? 2 + 0(g), the relation j^ 2 '(t) = 2j( 1 ' is not satisfied due to 
the e~ 2r+t factor and j ( 2 ) is exponentially smaller than 2j^' at t > 0. 

Let us compare the result with defined by eqn 
component of the local density operator is defined as 



P[0] (*) 
P[i] (t) 

P[2](t) 
P[l,2](*) 



r® 2 (g,t)[p(0)]=pf\t), 
T(g,t)[Tr 1 [p(0)}]-Tr 1 ' 
r(^,t)[Tr 2 [p(0)]]-Tr 2 



T^(g,t)[p(0)] 
T^(g,t)[p(0)] 



In this case each 

(129) 
(130) 
(131) 



1 - Tr 2 \r(g,t)[Tri\p(0)]]] - Tr x [T(g, t)[Tr 2 [p(0)}}} 



+Tr li2 T^(g,t)[p(0)] 



(132) 



As seen above, to calculate p^], p[ 2 ] and pnph we need the time evolution of 
the partial density matrices F(g, t)[Tri[p(0))) and T(g, t)[Tr 2 [p(0)]]. It should 
be stressed that these evolutions cannot be obtained only from knowledge of 
p^\t), for instance, by taking any traces forp^\ They must be calculated 
independently by solving eqn (jllfi|) for the initial density matrices Tri[p(0)] 
and Tr 2 [p(0)]. For the initial i.i.d. density matrix, each component is 
calculated as 



(133) 
(134) 
(135) 
(136) 



t(2) 
J M 


= j(2) = 8d(t) 2 e- 2r ^ 


t + 0(g) 


t(2) 
J [l] 


= Ad(t) 2 (l-e- 2r+t ] 


+ 0(g), 


t(2) 

J m 


= 4d(t) 2 (l - e- 2 ^] 


+ 0(g), 


7(2) 

J [l,2] 


= 0(g). 





Thus total information is precisely equal to twice of JW; 

JCO = 8d(t) 2 + 0(g) = 2JW. 
Next let us discuss an entangled case. Initially we take a state as 

l*(o)> = ^[| + -> + !-+>]. 

Calculation of j^(t) is easy and the results are as follows. 



(137) 



(138) 



3 W(t) = 8d(t) 2 [ 



2,-2F+t 



-2F-U 



+ 0(g). 



(139) 
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Note that j^(t)/j^(t) vanishes exponentially in time just as in the i.i.d. 
case. 

Evaluation of (t) needs not only the density matrix 

Pll (t) = i$(o> m)\ (wo) 

but also another density matrix 



T(g,t)[TrM0)]] = r(<M)[lY 2 [p(0)]] = F(g,t) 
After some manipulations the form of (t) results in 



1 

2 1 



;i4i) 



8rf(t) 2 [e- 2r +' + e- 2r -*] 
+4rf(t) 2 [l + 2e^ r++r -) t ] 2 — 



r_ r+t _ e -r_ t] 2 



e -zr+t(i _ e -2r_t) + e -2r.t^ _ e -2r +< ) 
+0(<?). (142) 

Note that at the early era (t ~ 0), both j^(t) and J^ 2 \t) have four-times 
information compared with the single system: 

j (2) (t~0) ~4j (1) , (143) 
J( 2 )(t~0) ~4J«. (144) 

Thus the information is twice larger than the above i.i.d. case. Obviously 
this advantage arises due to the entanglement between subsystems. 

For the entangled case, J^/jW becomes smaller than the value of the 
i.i.d. case (equal to two ) in the late time. Hence, the i.i.d. density operator 
becomes more relevant than the entangled density operator for the estimation 
of g. In the limit of t — > oo, the value of J^/J^ for the entangled case 
approaches to the single-system value: 

lim ^TTTT^ ~ !■ (145) 
t-ooj(l)(t) V 1 

This is due to contributions of the one-blank states (|-B±) and | ±B)). Con- 
sequently, it can be said that the measurement should be at the early times in 
order to utilize enhancement of the Fisher information by the entanglement. 
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So far we have discussed only systems with small numbers of samples. 
For the practical estimation of the small parameter g, the many-samples 
estimation is inevitable beyond the above simple examples. For instance, the 
minimized expected error 5g is given by 

5g = - 1 (146) 

for the i.i.d. cases of iV-samples systems. Then, in order to get a meaningful 
estimate, the number of the samples must be, at least, O (l/(g 2 J^)^j for 
the correct value g. It is expected that large entanglement between many 
samples may extremely improve the estimation for g and make the number 
of the samples enough for the estimation much smaller. 



11 Summary 



We have investigated deeply the local quantum estimation problem of an 
unknown parameter. The practical restriction of experimental observables 
takes place in various situations of the physical experiments. For a typical 
example, in particle physics we can probe only low-energy visible sectors of 
the whole system by our present devices. Such obstacles appear because 
of the limit of the present technology and so on. Moreover, observation of 
quantum phenomena, which happen only at quite small rates, often becomes 
the crucial target of experiments, which may derive some profound results of 
physics like, for instance, CPT violation [H]. In such situations, it is generally 
difficult to take a large number of data as one wishes, at least, in the first stage 
of the experimental studies. Hence the local quantum estimation becomes 
really important when the experimental arrangements are designed, because 
the estimation theory provides among our available probes the optimized 
observable which quantum fluctuation is most suppressed in the estimation 
based upon a limited number of the data. 

In spite of such relevance of the local quantum estimation, the problem 
has never been discussed in detail, as far as the authors know. In this paper, 
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the detailed analysis and formulations based upon the Fisher informations 
have been completed. After a brief review on the standard quantum estima- 
tion theory, the local quantum estimator for the local estimation has been 
defined by eqn (|42|) . The notion of the local density operators was clearly 
introduced in eqn ()47|) . and the Cramer- Rao inequality in the local quantum 
estimate theory (eqn (|61j) ) has been proven by taking the local Fisher infor- 
mation defined by eqn (J59j) . The inequality is a fundamental tool in the theory 
and will play a significant role in the local estimation in various physical appli- 
cations. In section 6, the Fisher information for the unnormalized pure state 
was commented. The formula in eqn()67|) is an extension of that derived by 
Fujiwara and Nagaoka, who discussed the Fisher information for normalized 
pure states. It is known that in many physical systems non-unitary theories 
of pure states also are available and that the validity is well verified by the 
experiments. In such systems with non-unitary evolution, eqn ()67|) is quite 
useful to evaluate the Fisher information for an unknown parameter. In sec- 
tion 7, it was pointed out that the local quantum estimation in the composite 
system has two independent formulations, using the i.i.d. cases. In section 8, 
two general formulations of the local quantum estimation for the composite 
system were proposed. For the composite system of N identical subsystems, 
we have two Fisher informations, and j( N \ The information j( N ) takes 
a simple form to define, but gives, in general, much smaller values than j( N '. 
The theory of information j( N ' can generate a more precise estimate for g, 
but has a pretty complicated form to deal with, compared to the j( N ' case. In 
order to avoid the troublesome procedures in evaluation of J( N > , we showed in 
section 9 the formula in eqn (jll2j) . which makes the evaluation tractable. As 
seen in eqn flll2|) . calculation of the Fisher information j( N > requires solving 



evolutions of many descendant operators, F^ N „ t £\(g,t) T^(ji,---,j m )\ptat 

independently of solving the accessible density operator £>y (t) itself. Such 
processes never appear in the ordinary quantum estimation theory, where the 
Fisher information can be evaluated by using only a time-evolved density op- 
erator. In section 10, we demonstrated explicitly our formulation of the local 
quantum estimation by applying to a decaying two-level system with a small 
unknown parameter. 

We hope that the analysis in this paper enables the quantum estimation 
theory to take a more active part in the real experimental studies, which suffer 
from the restriction of available observables and the practical limitation of 
the number of the data. 
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Appendix 

In this appendix, the Cramer-Rao inequality is proved. Let us write the 
triangular inequality relation as 

Tr(XiX)Tr(Y*Y) > \Tr(X^Y)\ 2 , (147) 

where X and Y are arbitrary operators acting on the Hilbert space. De- 
composing the operator X*Y into the sum of the real and imaginary parts 
as 

x f r = -(x + y + Y + x) + -{x^y-y^x), (148) 

2 2 
another inequality relation arises: 

Tr(X j X)Tr(Y j Y) > i|Tr(X f r + Y ] X)\ 2 + i|Tr(X f r - Y ] X)\ 2 

> ^\Tr(X^Y + Y t X)| 2 . (149) 

Here let us take 



X = L(g)J Ptot (g), (150) 



Y = (A-E g [A})^p tot (g). (151) 
Then, from the inequality (114911 . we can derive that 



TT[p tot (9)L(9) 2 mPtot(g)(A - E g [A]) 2 ] 

> \ |Tr \p tot {g) (L( 9 )(A - E g [A}) + {A - E g [A})L{g))}\ 2 . (152) 

The right-hand-side term in the above inequality is able to be calculated 
using eqns((TUJ), Q and (j3J) successively as follows. 

\ |Tr \p tot {g) (L(g)(A - E g [A)) + (A - E g [A])L(g))]\ 2 
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= \\Tr[p tot (g)(L(g)A + AL(g))}\ 2 

= ±\Tr[A( Ptot (g)L(g) + L(g)ptot(g))]\ 2 

= (Tr[Ad gPtat (g)]) a = (d g E g [A]f . (153) 
Consequently the relation ()152j) implies the following inequality: 

V ' lA] >|, (154) 



thus, the inequality (f2l)|) is proved. For the unbiased case with E g [A] = g, the 
inequality ()154j) is reduced to (pH?|) . The equality is trivially attained when 
X oc Y in eqn ()149p and the relation X oc Y holds in eqns ()150|) and (|151|) 
when we set A oc £(g), because E g [L(g)] = 0. 
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